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Convergence of density-functional supercell calculations for defect formation energies, charge tran- 
sition levels, localized defect state properties, and defect atomic structure and relaxation is inves- 
tigated using the arsenic split interstitial in GaAs as an example. Supercells containing up to 217 
atoms and a variety of k-space sampling schemes are considered. It is shown that a good description 
of the localized defect state dispersion and charge state transition levels requires at least a 217-atom 
supercell, although the defect structure and atomic relaxations can be well converged in a 65-atom 
cell. Formation energies are calculated for the As split interstitial, Ga vacancy, and As antisite 
defects in GaAs, taking into account the dependence upon chemical potential and Fermi energy. 
It is found that equilibrium concentrations of As interstitials will be much lower than equilibrium 
concentrations of As antisites in As-rich, n-type or semi-insulating GaAs. 

PACS numbers: 61.72.Bb, 61.72.Ji, 71.55.-i 
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I. INTRODUCTION 

Interstitials are the most complicated of the simple 
point defects, and the most elusive. For example, even 
though arsenic interstitials must be created by irradiation 
of GaAs with sufficiently energetic particles, and they can 
subsequently be observed to recombine with arsenic va- 
cancies when the sample is heated above 220° C, isolated 
arsenic interstitials have not been observed directly in 
EPR, electrical, or optical experiments^ 

It has been argued based on a thorough analysis 2 - - — 
of a variety of experimental data including titration 
experiments^ and measurements of density and lattice 
parameter— that melt-grown GaAs is always As-rich un- 
less the concentration of Ga in the melt is substantially 
greater than 50%, and that this deviation from stoichiom- 
etry is due primarily to the creation of large numbers of 
As interstitials (Asj) during growth. In particular, Hurlc 
has argued 2 - that the measured deviation of the mass per- 
unit cell as a function of arsenic concentration in the melt 
must be explained by arsenic interstitials and/or arsenic 
vacancies, since the number of arsenic antisites which 
would be required to fit the data is unrealistically large 
(up to several percent), due to the small difference be- 
tween the atomic masses of arsenic and gallium. Hurle's 
work also contains an extensive thermodynamic analysis, 
including estimates of the mass action constants for the 
formation of all the neutral native point defects. These 
estimates are derived by fitting to a large quantity of ex- 
perimental data on both doped and undoped GaAs, un- 
der the assumption that native defect and dopant concen- 
trations are near equilibrium close to the melting point 
and during high temperature growth from the melt or 
from solution. 

In the high temperature growth regime, observations 
of defects tentatively described as high concentrations 



or diffuse 'clouds' of arsenic interstitials have been re- 
ported in GaAs grown by the horizontal Bridgman and 
liquid-encapsulated Czochralski methods, based on X-ray 
diffuse scattering^— and quasi-forbidden X-ray reflection 
intensity measurements However, the atomic compo- 
sition and microscopic structure of these defects cannot 
be unambiguously determined from these experiments. 

Gallium arsenide grown by arsenic-rich molecular 
beam epitaxy at low temperature (LT GaAs) is a semi- 
insulating material with a host of potentially useful 
applications. 11 This material contains up to 1.5% excess 
Asji2 which is accommodated by high concentrations of 
point defects in UN-annealed samples, and arsenic pre- 
cipitates plus somewhat lower concentrations of point de- 
fects in annealed samples. Concentrations of As antisites 
(Asca) up to 10 20 cm' 3 are observed in LT GaAs, as 
measured by electron paramagnetic resonance (EPR), 13 
near-infrared absorption (NIRA) and magnetic circular 
dichroism of absorption (MCDA), 14 and scanning tunnel- 
ing microscopy (STM). 15 Concentrations of Ga vacancies 
(Vba) up to 10 18 cm -3 are measured in LT GaAs by slow 
positron annihilation; 1 - Ion channeling experiments have 
been interpreted as providing evidence for large concen- 
trations of As interstitials in LT GaAs. 12 However, it was 
later pointed out that the observed high concentration of 
atoms in the channel near the normal arsenic lattice sites 
could also be due to outward relaxation of the nearest 
neighbors of the As antisitesi 14 ' 17 

Within certain well defined limits of the growth param- 
eters for LT GaAs, a linear correlation between the neu- 
tral Asca concentration and the lattice dilation has been 
found. 14,16 It was therefore proposed that Asca are the 
dominant defects which determine the lattice expansion 
for growth within this regime. Staab et alM. used a self- 
consistent density-functional-based tight-binding method 
to study the lattice distortion induced by point defects 
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in As-rich GaAs, and concluded that only Asc a are nec- 
essary to understand the observed lattice expansion in 
the regime where the linear correlation is observed, and 
that if concentrations of isolated As^ comparable to the 
measured concentrations of Asc a were also present, the 
lattice expansion would be three times greater than is 
experimentally observed. However, Luysberg et al. have 
reported that when the As/Ga flux ratio is increased be- 
yond a beam equivalent pressure (BEP) ratio of 20, there 
is a departure from the linear correlation between lattice 
dilation and antisite concentration^ It was pointed out 
that other defects must be present to account for the de- 
viation from stoichiometry and the lattice expansion at 
high As/Ga flux ratios*^ 

Nonequilibrium processes such as diffusion and compo- 
sitional intermixing at interfaces can also be strongly af- 
fected by point defects that are present in high concentra- 
tions. Since the point defects which have been unambigu- 
ously documented as present in high concentrations in LT 
GaAs, Asca and Vca, occupy sites on the gallium sub- 
lattice, they cannot contribute directly to intcrdiffusion 
on the arsenic sublattice. However, substantial concen- 
trations of arsenic intcrstitials may affect intcrdiffusion 
on the arsenic sublattice. For example, an experimental 
study showing a positive dependence of GaAsP/GaAs 
and GaAsSb/GaAs interdiffusion on arsenic pressure has 
indicated that a kickout mechanism involving arsenic in- 
terstitials is the dominant process for the As-P and As-Sb 
intcrdiffusion in the material studied*^ 

Similarly, annealing LT-GaAs 5-doped with Sb was 
found to produce substantially greater compositional 
intermixing than annealing conventional stoichiometric 
GaAs similarly (5-doped; 1 - This enhancement of As-Sb 
interdiffusion was attributed to an oversaturation of ar- 
senic interstitials in the LT GaAs samples, resulting from 
the balance of arsenic interstitials with arsenic clusters 
and all the other excess-arsenic-containing defects in the 
material. The effective activation energy for As-Sb inter- 
diffusion in LT GaAs deduced from this work, 0.6 ± 0.15 
eVji£ is reasonably close to the migration energy of 0.5 
eV for arsenic interstitials deduced from annealing ex- 
periments on defects produced by electron irradiation, 1 
as well as to migration energies subsequently ascribed 
to arsenic interstitial defects produced in GaAs by other 
means. The concentration of arsenic interstitials required 
to produce sufficient oversaturation to eliminate com- 
pletely any contribution of the interstitial formation en- 
ergy to the activation energy for As-Sb intermixing mea- 
sured in the LT GaAs sample was estimated to be roughly 
10 18 cm -3 , using Hurle's thermodynamic analysis in con- 
junction with the experimental dataJ^ 

Theoretical attempts to obtain a picture of the micro- 
scopic structure and properties of the lowest energy ar- 
senic interstitial configuration(s) began with the work of 
Baraff and Schliiter, who used density functional Green's 
function calculations to investigate the energies of re- 
actions creating native defects with Td symmetry in 
GaAs, including arsenic interstitials in the two tetrahc- 



dral sites^ The effects of lattice relaxation were ignored. 
Baraff and Schliiter concluded that simple tetrahedral ar- 
senic interstitials were less likely to occur than vacancy 
and antisite defects under all equilibrium conditions, al- 
though they could not rule out the possibility that other, 
more complicated interstitial configurations might have 
a lower energy^ 

Jansen and Sankey calculated the formation energies 
for unrelaxed native defects with tetrahedral symme- 
try in GaAs, including arsenic interstitials in tetrahedral 
sites, using a density-functional pseudopotential method 
with a basis set of pseudo-atomic orbitals and a single 
special k-point in supercells containing about 32 atoms*2A 
In order to calculate the formation energies for individual 
defects instead of reaction energies for defect reactions 
which conserve the number of atoms of each species, they 
were required to choose a value for the arsenic chemical 
potential (or equivalently, for the gallium chemical po- 
tential) . An arbitrary value was chosen, corresponding to 
the condition that the formation energies for neutral gal- 
lium vacancies and for neutral arsenic vacancies should 
be equal. Jansen and Sankey concluded that arsenic in- 
terstitials in tetrahedral sites should be less numerous 
than vacancies and antisites in GaAs under equilibrium 
conditions,-— in agreement with Baraff and Schliiter. 

Zhang and Northrup used density functional theory 
(DFT) within the local density approximation (LDA) 
and supercells of about 32 atoms to calculate the for- 
mation energies for vacancies, antisites, and tetrahedral 
interstitials in GaAs as a function of arsenic chemical 
potential, over the physically allowable range of chemi- 
cal potentials, from Ga-rich to As-rich^ This physically 
allowable range is set by the heat of formation of bulk 
GaAs and by the requirement that the arsenic chemical 
potential may not exceed the chemical potential of bulk 
arsenic, since the material is in equilibrium with arsenic 
precipitates in the arsenic-rich limit. The atomic coordi- 
nates were allowed to relax in these calculations, within 
the constraints imposed by the tetrahedral symmetry. In 
agreement with the previous work, Zhang and Northrup 
found that antisites and/or vacancies should be more nu- 
merous than arsenic interstitials in tetrahedral sites un- 
der all equilibrium conditions^ 

Chadi used DFT-LDA calculations and 33-atom super- 
cells to investigate many different types of bonding con- 
figurations for sclf-intcrstitials in GaAs, including vari- 
ous split interstitials, as well as hexagonal, two-fold co- 
ordinated, and tetrahedral interstitials, all fully relaxed 
within the constraints of the chosen symmetry^ He 
found that the lowest energy configuration for arsenic 
interstitials in the neutral or —1 charge state is a split 
interstitial consisting of two As atoms sharing an arsenic 
lattice site, displaced from this site in opposite directions 
along a <110>-like axis, while the lowest energy config- 
uration for positively charged arsenic interstitials in the 
+1 or +2 charge state is a split interstitial consisting of 
an As atom and a Ga atom sharing an gallium lattice 
site, displaced from this site in opposite directions along 
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a <100>-like axis. Since we will be interested below pri- 
marily in arsenic interstitials in semi-insulating or n-type 
GaAs, we will use the notation As^-As for the interstitial 
with two atoms sharing an arsenic site and aligned along 
a <110>-like axis, which should be the lowest energy in- 
terstitial configuration in semi-insulating or n-type ma- 
terial. 

Chadi also showed that neutral arsenic interstitials, 
which have unpaired spins, are unstable relative to for- 
mation of a pair of +1 and —1 charged interstitials - 
i.e. arsenic interstitials form a negative-f7 system. This 
suggested that arsenic interstitials may not be observ- 
able in EPR experiments* 2 ^. Chadi reported the relative 
energies for the most energetically favorable arsenic in- 
terstitial configurations in each of these charge states, 
including in each case a number of metastable configura- 
tions somewhat higher in energy than the lowest energy 
configurations, all of which were more complicated than 
the simple tetrahedral configurations. 23 However, since 
Chadi did not report absolute interstitial formation en- 
ergies as a function of arsenic chemical potential, no com- 
parison with the formation energies of defects involving a 
different number of excess arsenic atoms, such as arsenic 
antisites, was possible from this work. 

Landman et al. investigated the relative formation en- 
ergies of the point defects containing excess As, Asca, 
Vca, and the lowest energy As^ configuration in semi- 
insulating or n-type, As-rich GaAs, As^-As* 2 ^ using DFT- 
LDA-based calculations with the Harris-Foulkes func- 
tional and a basis of pseudo-atomic orbitalsj 25 ' 26 They 
placed the defects in 64-atom supercells, and estimated 
summations in k-space by using a single Chadi and Co- 
hen special point, 27 Since Harris-Foulkes, pseudoatomic- 
orbital calculations do not give as accurate results for 
semiconductor heats of formation or for the relative 
energies of compound semiconductor and pure, metal- 
lic phases as fully self-consistent DFT-LDA calculations 
with a sufficiently large basis of plane waves, they did not 
use this method to calculate the arsenic chemical poten- 
tial in the arsenic- rich limit. Instead, the relative forma- 
tion energies of the tetrahedral As^ with arsenic nearest 
neighbors, Asca, and Vba m the arsenic-rich limit for the 
chemical potential were taken from information given by 
Zhang and Northrupj 2 ^ and the relative formation energy 
of As^-As was determined by Landman et aL's result that 
the neutral As^-As is 4 eV lower in energy than the unre- 
laxed, neutral tetrahedral As^ with arsenic nearest neigh- 
bors. Since the tetrahedral interstitial was found to be 
unstable, relaxing to another configuration in Landman 
et aZ.'s calculation, they were obliged to compare their 
results for the ideal, unrelaxed tetrahedral interstitial to 
Zhang and Northrup's results for a tetrahedral intersti- 
tial which had been relaxed while constrained to keep 
its tetrahedral symmetry. This led to an additional un- 
certainty in the relative formation energies between zero 
and 0.8 eV<2i However, Landman et al. concluded that 
the lowest energy split As^ may have a concentration ap- 
proaching that of Asca for certain Fermi levels^ 



Since the theoretical investigations described above 
have been carried out over a long period of time, it has 
gradually become possible not only to include lattice re- 
laxation and to investigate more complicated interstitial 
configurations, but also to do more accurate calculations, 
using larger unit cells and better sets of k-points for the 
summations over k-space. Poykko et al. showed how sen- 
sitive calculated defect properties can be to the k-space 
sampling method and supercell size in their investigation 
of the VAs-SiGa complex in GaAs^ They found that the 
use of the T point can produce misleading results even 
when supercells are 64-atoms in size, reinforcing the con- 
clusions of Makov that the T point produces particularly 
slowly converging results with respect to cell size^ So it 
is essential to use a special point mesh in this type of cal- 
culation. Furthermore, Puska et al. concluded that cell 
sizes of 128 to 216 atoms are needed to properly assess 
the physical properties of the silicon vacancy in bulk sil- 
icon, because of the dispersion of energy levels and long 
range ionic relationships^ 

In this paper, we investigate the combined effects of 
cell size and k-space sampling on the formation energy, 
charge state transitions, atomic relaxations, and char- 
acterization of localized defect states for arsenic self- 
interstitials in GaAs. Because of the more ionic nature of 
the material and the complicated split interstitial defect 
structure, comparison of these results for interstitials in 
GaAs to the previous results for vacancies in silicon 3 -^ can 
enhance our understanding of the range of behavior for 
different defects in different materials. We compare the 
formation energy of the lowest energy arsenic interstitial 
in n-type or semi-insulating GaAs, Asj-As, with the for- 
mation energies of AsQ a and Vqh at the arsenic-rich end 
of the range of physically allowed chemical potentials, 
all calculated by state-of-the-art DFT pseudopotential 3 ^ 
calculations, using the larger supercells and sets of spe- 
cial k-points which we have determined to be necessary. 
We conclude our study by discussing the relative con- 
centrations of these defects in equilibrium in As-rich, n- 
type or semi-insulating GaAs at growth temperatures, 
and reporting the computed charge transition levels and 
expected electrical behavior of As,-As as a function of 
Fermi level. 



II. COMPUTATIONAL METHOD 

We have used the molecular dynamics code developed 
at the Fritz Haber Institut (FHIMD) 31 for this investiga- 
tion, using density-functional theory (DFT) 32 within the 
local density approximation (LDA), with the Ceperley- 
Aldcr 3 — form for the exchange and correlation poten- 
tials as parameterized by Perdew and Zunger^ The 
core electrons are treated in the frozen-core approxima- 
tion and the ion cores are replaced by fully-separable 3 ^ 
norm-conserving pseudopotentials. 36 Plane waves are in- 
cluded up to the energy cutoff of 10 Ry. The atoms 
are allowed to relax until the force components are are 
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less than 5 x 10 ~ 4 hartrees per bohr radius and the 
zero temperature formation energies change by less than 
5 x 10~ 6 hartrees per step for at least 100 steps. 

To evaluate the defect formation energy, we used the 
formalism of Zhang and Northrup, 22 which gives for the 
formation energy in the As-rich limit at zero temperature 

AE f = E{N Gai N Asi q) - A Ga /i GaAs 

-(^As - A Ga )/i As(bulk) + Q^F ■ (1) 

Here q electrons have been transferred to a reservoir at 
the Fermi energy cf in order to produce a defect in 
the desired charge state. E(Ngh, Nas> q) is the zero- 
temperature total energy produced by the ab initio code 
for a supercell containing the desired defect, the chemi- 
cal potential /XQaAs is the energy per atomic pair of bulk 
GaAs, and the arsenic chemical potential in the As- rich 
limit, /XAs(buik)) is the energy per atom of pure bulk As 
computed using the same ab initio code and pseudopo- 
tentials. Nq & and Na s are the numbers of atoms of each 
species in the supercell containing the defect. We will dis- 
cuss the effect of temperature, which can be important 
for defect concentrations, in Section lull 

Because the zero of the energy levels floats freely^ re- 
sults from different DFT supercell calculations must be 
aligned in order to obtain the correct charge transition 
levels. We apply the procedure outlined by Kohan et 
al£^ in which we first compute the difference between 
the electrical potential in the supercell with the neutral 
defect and the electrical potential for the corresponding 
bulk crystal supercell, averaged over parallel planes, as 
a function of position along a line normal to the planes. 
Far from the defect within the supercell, this difference 
becomes a constant. In order to make the potential far 
from the neutral defect equal to the corresponding po- 
tential in the bulk cell, a uniform shift is applied to the 
potential and the energy levels, yielding the proper align- 
ment of the energy levels of the defect with the energy 
levels of ideal bulk crystal supercell. The same shift is 
applied for all charge states of a given defect. 

A well-known shortcoming of the LDA is that it under- 
estimates the band gaps of materials. The typical method 
for dealing with this problem is to simply shift the con- 
duction band states up uniformly by the amount needed 
to reproduce the experimental gap, using the so-called 
'scissors operator'.— More recently, an analytically-based 
model justifying the rigid shifting upward of all conduc- 
tion band states by a scissors-type correction has been 
shown to produce good results for a large number of 
semiconductors.— Since the LDA can produce similarly 
large errors in the energies of the deep defect states, it 
is also important to correct for these errors when deter- 
mining where the charge transition levels corresponding 
to deep defect states lie in the experimental gap. Unfor- 
tunately, a full GW calculation^ - — which would correct 
these errors, is not currently possible for the large super- 
cells needed for studies of defects. Therefore, we apply 
the same upward shift to the defect states with predom- 
inantly conduction band character as is applied to the 



conduction band states themselves, while leaving the de- 
fect states with predominantly valence band character 
fixed relative to the valence band edge. 



III. RESULTS AND DISCUSSION 

A. Defect formation energies, charge transition 
levels, and localized defect states 

We used cubic supercells with dimensions of both two 
and three times the computationally determined bulk lat- 
tice constant, corresponding to bulk cells of 64 and 216 
atoms, along with three different Brillouin zone (BZ) 
sampling schemes to examine the effects cell size and 
sampling scheme have on the formation energies and 
transition levels for the As^-As. In order to investigate 
the most efficient choice of k-points to obtain good re- 
sults, we have used a l 3 Monkhorst-Pack (MP) mesh^ a 
2 3 MP mesh, and the T and L points, which was recom- 
mended as a good minimal set of k-points for cubic super- 
cells with no particular symmetry within the supercell ^ 
and which has subsequently been used in defect calcula- 
tions, e.g. to study the structures associated with dopants 
in highly n-doped Si.~ Calculations comparing different 
summation schemes for vacancies in Si show that use of 
the r + L points produces a reasonably well-converged 
formation energy in a 64-atom supercell. 30 

This earlier work on the vacancy in silicon has shown 
that the neutral vacancy formation energy computed 
with different k-space sampling schemes converges at dif- 
ferent rates with respect to supercell size4& However, an 
acceptably converged value for the neutral vacancy for- 
mation energy can be attained more easily than an ac- 
ceptably converged description of the charge transition 
levels and the atomic relaxations and defect symmetry 
for different charge states.— 

To augment our understanding of the effects of cell 
size and sampling scheme, we computed the formation 
energies of the different charge states for the fully relaxed 
split interstitial As^-As in GaAs. Charge was balanced 
by a uniform background to avoid long range Coulomb 
interactions between the supercells. 

Table Q] lists the formation energies we obtained for all 
the excess-arsenic-containing elementary point defects in 
the arsenic-rich limit, with GaAs in equilibrium with bulk 
arsenic, including complete results for the various k-space 
sums for Asi-As. For comparison, the formation energies 
of the unrelaxed, ideal tetrahedral As interstitials with 
As neighbors (Asji) and with Ga neighbors (Asi2) are 
also shown. These tetrahedral interstitials are unstable, 
and will relax to other configurations if allowed to break 
their tetrahedral symmetry. 

Table |T] displays formation energies evaluated for the 
Fermi level at the valence band maximum (VBM), and 
also for the Fermi level pinned at the calculated position 
of the (+1/0) charge state transition of the Asca, which is 
at VBM+0.54 eV for the 2 3 MP mesh and VBM+0.45 eV 
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TABLE I. Formation energies for excess-arsenic-containing 
defects computed in the As-rich limit. These were calculated 
with supercells corresponding to the bulk 216 atom cell. The 
values in the last column are computed with the Fermi level 
pinned at the calculated (+1/0) transition level of the As an- 
tisite defect. 
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FIG. 1. Defect formation energies for the As split interstitial 
computed in the As-rich limit. Comparisons are presented 
for different k-space sums and two supercell sizes. Dashed 
lines are used for cells containing 65 atoms, and solid lines for 
217 atom cells. Transition levels between the different charge 
states are marked with circles for the 2 3 MP mesh, squares 
for the l 3 MP mesh, and triangles for the T + L k-space sum. 



for the l 3 MP mesh. This choice of Fermi level was based 
on the experimental finding that there can be high con- 
centrations of AsQa m As-rich GaAs, and that these high 
concentrations of Asc a can pin the Fermi energy near 
this transition level. Using either choice of Fermi level as 
a reference, the formation energies and equilibrium con- 
centrations of the defects can be determined as a function 
of Fermi energy (or doping level) . 

In Fig. [TJ we present the results of the formation en- 
ergy calculations for Asj-As in which both cell size and 
k-space sampling methods have been varied. This figure 
shows the formation energy for the Asi-As in its pre- 
ferred charge state. For Fermi levels in the range where 
the neutral As^-As is preferred, the formation energy is 
independent of Fermi level. If the Fermi level is decreased 
past the (0/+1) charge transition level, so that the +1 
charge state is preferred, the formation energy vs. Fermi 
level curve has a slope of +1, as can be seen from Eq. ((T|), 
and if the Fermi level is increased past the (0/— 1) charge 
transition level, so that the —1 charge state is preferred, 
the formation energy vs. Fermi level curve has a slope of 
-1. 

The most obvious feature in Fig. [1] is the wider varia- 
tion in the energies computed using the smaller supercell. 
Convergence with cell size is also visibly slower when us- 
ing the l 3 MP mesh or the T + L points, which use a 



less finely spaced set of k-points than the 2 3 MP mesh to 
cover the Brillouin zone. It is clear that the 217 atom cell 
with the 2 3 MP mesh is well converged. We see that the 
two less finely spaced sampling schemes produce some- 
what converged results in the 217 atom cell, as does the 
2 3 MP mesh in the 65 atom cell. In agreement with pre- 
vious results for the vacancy in bulk Sij^H we find that 
different sampling schemes can be used to produce either 
an attraction between defects (i.e. a lowering of energy 
in the smaller supercell), as in the l 3 MP or the T + L 
cases, or a charge-state dependent repulsion or attraction 
between defects, as in the 2 3 MP case. 

When comparing the results for the 65-atom and the 
217-atom cells, using the 2 3 MP mesh, we see that the 
charge transition levels converge less rapidly than the 
neutral defect formation energy with increasing cell size, 
due to long-range interactions of the electrons in the lo- 
calized defect state with the charged defects in neighbor- 
ing cells. The neutral As^-As formation energy changes 
by less than 0.1 eV when the cell size is reduced from 217 
to 65 atoms, while the charge transition levels move by 
about 0.4 eV, when using the 2 3 MP mesh. This trend is 
not seen in the less- well-converged results obtained using 
the l 3 MP mesh or the T + L points, where the neutral 
As^As formation energies move by 0.5 to 0.6 eV when 
the cell size is reduced from 217 to 65 atoms, perhaps 
due to interactions of the deep defect level with the band 
edges. 
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FIG. 2. This contour plot shows levels of constant charge 
density for the localized defect state of the neutral Asi-As 
in the 217 atom supercell, evaluated in the plane containing 
the two defect As atoms, which is 0.48A above the As lattice 
site associated with the defect. The dark circles represent the 
positions of the As atoms, including the defect atoms and the 
As atoms of the original lattice plane. The light circles show 
the locations of the neighboring Ga atoms projected onto the 
plane. 



In the neutral charge state, the topmost filled elec- 
tronic level is half filled. This level corresponds to a local- 
ized state centered on the split interstitial oriented along 
a <110>-like direction. In Fig. [21 the charge density as- 
sociated with this defect state is shown in a plane parallel 
to an arsenic lattice plane, but slightly above this plane, 
so that it includes the two arsenic atoms of the split in- 
terstitial, which have relaxed slightly away from the ideal 
lattice site. This plot clearly shows that the defect state 
is localized in p-like orbitals which appear to be forming 
a 7r antibonding state, as evident from the node in the 
charge density midway between the two As atoms form- 
ing the split interstitial. This result is corroborated by 
an examination of the characterization of the state, using 
the 217-atom cell and the 2 3 MP mesh, which shows that 
in contrast to the very extended character of the bulk- 
like states, 20% of the defect state is localized in the p y 
and p z orbitals of the two As atoms of the neutral de- 
fect. Since the charge density in this defect state points 
roughly along the direction of the As-Ga bonds between 
the defect As atoms and the two Ga atoms bonded to 
both of these As atoms, we note that this state makes a 
bonding contribution to these As-Ga interactions as well 
as making an antibonding contribution to the interaction 
between the As atoms of the defect. 

The characterization of the deep defect state in the 
217-atom supercell is not very sensitive to the k-space 
sampling method used. For example, 20% of the state 



is also found to be localized in the p y and p z orbitals of 
the two As atoms of the neutral defect when the T + L 
sampling is used. 

In smaller supercells, the defect state interacts with its 
images and becomes less localized. One way to observe 
this is to examine the variation in energy of the defect 
state across k-space, i.e. the dispersion of the state. We 
characterize this variation by computing the difference 
between the highest and lowest energy obtained from a 
k-space survey along the A, £, A, and T lines in the cubic 
Brillouin zone. The dispersion measured in this way is 
0.1 eV in the 217 atom supercell, and 0.5 eV in the 65 
atom cell. This interaction of the defects in neighboring 
supercells contributes to the variation in the positions 
of the charge transition levels between the different cell 
sizes and sampling schemes, noted above, and supports 
the conclusion that the larger 217 atom supercell should 
be used for accuracy in describing the charge transition 
levels and deep defect states. 



B. Defect atomic structure and relaxation 

The detailed structure of the atomic positions is also 
expected to be dependent on the cell size and k-space 
sampling approach. However, in contrast to the behav- 
ior observed for the vacancy in silicon^ where the defect 
symmetry and atomic relaxations are very sensitive to 
the supercell size and k-space sampling, we find that the 
atomic structure of the As^-As defect is remarkably simi- 
lar for the 65-atom and 217-atom supercells, and also de- 
pends little on the k-space sampling approach. In Fig. [3] 
we show the structure of the neutral split interstitial ori- 
ented along the <011> direction, viewed from the [100] 
direction. The two Ga atoms labeled Ga(l) and Ga(3) 
are bonded to both As atoms of the defect, while those 
labeled Ga(2) and Ga(4) are bonded to only one of the 
defect atoms. The structure exhibits Ci v symmetry in 
all charge states and returns to this symmetry when the 
atomic coordinates are perturbed from their equilibrium 
positions and allowed to relax. 

Although the local lattice expansion introduced by the 
additional atom of the interstitial might be expected to 
converge slowly with supercell size, we find that reason- 
able convergence is more easily reached for the atomic 
positions than for the electronic properties discussed in 
the previous section. The As-As and As-Ga distances 
found in the 65 atom supercell differ by less than 0.01 A 
from those found in the 217 atom supercell, using the 2 3 
MP mesh. The distance between the two As atoms of the 
Asi-As defect is 2.39 A. (For comparison, the GaAs bulk 
nearest neighbor distance in this calculation is 2.41 A.) 
The distance between Ga(l) or Ga(3) and either of the 
As atoms of the defect is 2.60 A. The distance between 
Ga(2) or Ga(4) and the nearest As of the split interstitial 
is 2.32 A. 

The two different supercell sizes produce slightly dif- 
ferent results in the position of the center of mass of 
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FIG. 3. The neutral As^-As defect from the 217-atom su- 
percell is shown (viewed from a direction slightly displaced 
from the [100] direction) with larger dark spheres represent- 
ing the As atoms and smaller light spheres representing the 
Ga atoms. 



the two As atoms of the split interstitial. The center of 
mass of the pair of As atoms is shifted away from the 
bulk lattice site by 0.48 A in the [100] direction toward 
the plane containing Ga(l) and Ga(3) in the case of the 
217 atom supercell, using the 2 3 MP mesh. This shift is 

0. 50 A in the smaller supercell. Since the nearest neigh- 
bor distances are very similar in the two supercells, this 
difference is accomplished through variation of bond an- 
gles. The angle between the bonds from one of the As 
atoms to the Ga(l) and Ga(3) atoms is 107.7° in the 
small supercell and 106.8° in the large supercell. 

We find that the atomic structure of the As^-As defect 
depends very little on the k-space summation method, 
as can be shown by examining the bond lengths between 
the As atoms of the defect and between those As atoms 
and the four neighboring Ga atoms in the 217 atom su- 
percell. For the neutral As^-As defect, these bond lengths 
change by less than 0.1% when the summation method 
is changed. 

For the +1 charge state of the defect, we find a small 
but perceptible dependence of the bond lengths on the k- 
space summation method. The As- As distance obtained 
using the 2 3 MP mesh is about 1% smaller than that 
obtained using the T + L sampling. The bond between 
either of the As atoms and the Ga(l) or Ga(3) atom 
is about 1% longer when using the 2 3 MP mesh than 
when using the T + L summation method. The distance 
between the Ga(2) or Ga(4) atom and the nearest As 
atom of the split interstitial is about 0.5% longer for the 
2 3 MP mesh calculation than for the T + L calculation. 

There is a similar but weak effect (under 0.5%) in the 
defect bond lengths observed in the —1 state, with the 
roles of the 2 3 MP mesh and the T + L points reversed — 

1. e. the 2 3 MP mesh now gives a larger As- As distance 
and smaller As-Ga distances. The l 3 MP mesh produces 
results between those of the other two k-space summation 
methods. 



The effect of these small variations in bond lengths is 
to reduce the changes in bond length seen in the 2 3 MP 
mesh calculations when the defect becomes charged, if 
one of the other k-space sampling methods is used in- 
stead. In particular, the T + L sampling is observed to 
produce a somewhat smaller dependence of these bond 
lengths on charge state. This dependence of bond lengths 
on the charge state of the defect is seen below to result 
from the changes in occupation of the deep defect state 
when the charge state is changed, and from the bond- 
ing or antibonding character of this state for particular 
bonds. 

Focusing on results for the 217 atom cell with the 
2 3 MP mesh, we observe that the distance between the 
atoms of the defect has a significant dependence on the 
charge state of the system. The As- As bond expands to 
2.47 A (about 3.5% compared to the neutral state bond 
length), when the system is allowed to relax in the — 1 
charge state. This can be easily understood, since the 
antibonding defect state on the two As atoms is doubly 
occupied and can cause a stronger repulsive contribution 
to the interaction between these two atoms for the — 1 
charge state of the defect, while it is only singly occupied 
for the neutral state of the defect. 

When the system is allowed to relax in the +1 state 
the As- As bond shrinks to 2.31 A, a contraction of about 
3%. This can also be easily understood, since the anti- 
bonding defect state on the two As atoms of the defect is 
completely unoccupied for the +1 charge state, so it no 
longer makes a repulsive contribution to the interaction 
between the two As atoms. The bond length between 
Ga(l) or Ga(3) and either As atom of the defect is 3.5% 
longer and the bond length between Ga(2) or Ga(4) and 
the nearest defect As atom is 1.9% longer in the +1 state 
than in the neutral state. Since the deep defect state acts 
as a bonding state for the As-Ga bonds between the de- 
fect As atoms and Ga(l) and Ga(3), as discussed above 
in Section UlI A[ it is easy to understand why these bonds 
are longer in the +1 state, when the defect state is fully 
unoccupied and can no longer contribute to the strength 
of these bonds. 

We note that because of the contribution of the defect 
state to the As-Ga bonds between the defect As atoms 
and Ga(l) and Ga(3), these two Ga atoms move signif- 
icantly when the charge state is changed, changing the 
occupation of the defect state. These Ga atoms are about 
4% closer to each other in the —1 state and about 5.5% 
farther apart in the +1 state than in the neutral state. 
The other two Ga atoms, each bonded to only one defect 
As atom, do not move in response to the change in charge 
state. 

In performing these calculations, we fixed the lattice 
constant at the value determined through minimization 
of the energy of the bulk crystal. While the ideal calcula- 
tion should include a full lattice constant determination 
with each change of defect configuration, for simplicity we 
did not perform this relaxation. This may be deemed a 
reasonable choice in light of evidence presented by Puska, 
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et al£Q for ab initio supercell calculations in the LDA, 
using supercells of sizes comparable to ours, in which va- 
cancies in bulk Si are found to alter the lattice constant 
by around 0.2%, while artificially introduced distortions 
in the lattice constant of up to 1% are seen not to affect 
their reported results significantly. 



C. Relative defect concentrations in equilibrium 

We now compare our well-converged results for the 
formation energies of the elementary excess-arsenic- 
containing point defects, Asc a , Vba, and As^-As (the 
most favorable As^ configuration in semi-insulating or n- 
type GaAs), computed using the large supercell and the 
2 3 MP mesh. These formation energies in the As-rich 
limit, corresponding to GaAs in equilibrium with bulk 
arsenic, are presented as a function of Fermi energy in 
Fig. |U The formation energies for two specific choices 
of Fermi energy have also been listed in Table U in Sec- 
tion HUH 

In Fig. 21 we can see that the Vca and Asca defects 
possess significantly lower formation energies than the 
As^-As for all Fermi energies. For example, the forma- 
tion energy for Asca is seen to be at least 2 eV lower 
than the formation energy for As^-As for all Fermi ener- 
gies. Small uncertainties in the formation energy should 
not alter this strong qualitative ordering of the formation 
energies or the prediction based on this ordering that 
equilibrium concentrations of As^-As should be signifi- 
cantly lower than equilibrium concentrations of Asca, as 
discussed below. 

To estimate equilibrium concentrations of the excess- 
arsenic-containing defects we begin with the usual ex- 
pression 

C = N e -^E f /k B T e Sf/k Be -PAV f /k B T ] , 2 ) 

where N is the number of sites at which the defect can 
form in the crystal per unit volume, AEf is the total 
energy of formation of the defect, ks is the Boltzmann 
constant, and T is the temperature. The formation en- 
tropy of the defect is Sf, P is the pressure, and AVf 
is the change in the crystal volume associated with the 
defect formation. 

We note that the defect formation energy AEf for 
charged defects, as given by the formula in Eq. (TTJ), has an 
explicit dependence on the Fermi energy, in addition to 
its dependence on the calculated energies for defect for- 
mation at zero Fermi energy. Therefore we must compute 
the Fermi energy self-consistently, in order to determine 
the native defect concentrations present in a particular 
sample. If any electrically active impurities or dopants 
are present in the material, the concentrations of these 
impurities or dopants in all charge states must also be 
taken into account. We must set up the charge balance 
equation, requiring that the free electron and hole con- 
centrations (which also depend on the Fermi level) must 
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FIG. 4. Defect formation energies for selected defects over 
the calculated band gap in the As-rich limit. Zero Fermi level 
corresponds to the valence band maximum. 

cancel out any net charge resulting from the concentra- 
tions of all positively and negatively charged defects and 
impurities. This equation can then be solved to deter- 
mine the Fermi level. Once the Fermi level is known, 
it may be used to determine the formation energies and 
the resulting equilibrium concentrations of all the defects 
present. 

If the defect formation energy Ef for the most ener- 
getically favorable charge state of As^-As in a particu- 
lar sample is within ksT of the formation energy of the 
most favorable charge state of Asca, we may expect that 
the equilibrium concentrations of these two defects are 
comparable, assuming that the effects of the entropy of 
formation Sf and the change in the crystal volume AVf 
associated with the defect formation can be neglected. 
We will now concentrate on the relative defect concentra- 
tions at 1500 K (near the melting point of GaAs), since 
defects with a higher formation energy such as Asi-As 
have their greatest chance to attain equilibrium concen- 
trations comparable to those of more energetically favor- 
able defects at high temperature. We will estimate the 
effective corrections to the formation energy which occur 
at this temperature due to the entropy of formation and 
change in volume associated with the defects. 

First, to estimate the effect of the change in volume, 
we let P be atmospheric pressure and overestimate AVf 
to be the volume per bulk atom in the cell, which gives an 
effective correction to the defect formation energy PAVf 
of 9 x 10 -5 eV. We may safely neglect this correction. 

Previous calculations^ on defects in Si found that the 
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formation entropy Sf is dominated by vibrational con- 
tributions, and that the formation entropies are and 
5A;_b for the self-interstitial and the vacancy, respectively. 
We note that the self-interstitial in silicon is a <110> 
split interstitial with the same basic structure as Asi-As 
in GaAs. Therefore, in analogy to the results for de- 
fects in silicon^ we may assume that it is unlikely for 
the split interstitial Asi-As to have a very different for- 
mation entropy when compared to defects such as Asc a , 
which only contain atoms occupying lattice sites. If we let 
Sj = 10ks (an overestimate) for Asj-As, this gives rise 
to an effective reduction of the defect formation energy 
by SfT, or 1.3 eV at 1500 K. Even if we apply no reduc- 
tion to the Asca formation energy due to entropy, this 
still leaves the effective formation energy about 0.7 eV 
higher for As^-As than for Asca, producing equilibrium 
concentrations of As^-As which are about 0.4% those of 
As Ga at 1500 K. 

We conclude that even using this extremely liberal es- 
timate for the formation entropy of As^-As and ignoring 
the formation entropy of AsQ a cannot lead to an As^-As 
concentration approaching that of the antisites in ther- 
mal equilibrium. 
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FIG. 5. Transition levels computed with a rigid shift applied 
to the conduction band and conduction band derived states 
to correct the LDA band gap underestimate. 



D. Defect electrical behavior 

Although the placement of the calculated charge tran- 
sition levels in the experimental gap has an uncertainty 
far exceeding 0.1 eV due to the shortcomings of the LDA 
in calculating the gap, as discussed in Section UH our as- 
calculated band gap and charge transition levels are re- 
ported here to 0.1 eV (or 0.01 eV, for the closely spaced 
Vcja levels) , for the convenience of the reader who prefers 
not to read them off the picture in Fig. 01 For the Asca, 
the (+2/+1) transition level appears at 0.4 eV above the 
VBM, and the (+1/0) level is at E VB m + 0.5 eV. For the 
Asi-As, the (+1/0) transition is at -EVbm + 0.3eV, and 
the (0/— 1) transition is at -EVbm + 0.5eV. The levels for 
the Vba defect are at 0.09 eV, 0.13 eV, and 0.2 eV above 
the VBM for the (0/-1), (-1/-2), and (-2/-3) transi- 
tions, respectively. The calculated band gap of 0.8 eV is 
underestimated by 0.7 eV compared to the experimental 
zero-temperature gap of 1.5 eV. 

As discussed previously in Section UH we can get a 
rough estimate of where the charge transition levels fall 
within the experimental gap by shifting the conduction 
band derived states (including the deep defect states with 
primarily conduction band character) by the amount 
needed to correct the gap, while leaving the defect states 
with predominantly valence band character fixed relative 
to the valence band edge. Since the acceptor levels of the 
Vca are derived from the dangling bonds on the arsenic 
neighbors of the vacancy, which require three extra elec- 
trons to fill them, they should have the predominantly 
valence band character of anion dangling bond states. In 
Section IIII A[ the deep defect state of the As^-As was 
shown to have predominantly arsenic p-type character, 



similar to the character of the valence band edge states. 
However, the AsQa double donor defect state derives from 
an antibonding state of predominantly conduction band 
character, which has been lowered in energy due to the re- 
placement of the original anion-cation bonds of the ideal 
crystal by anion-anion bonds between the antisite and its 
nearest neighbors. This donor state is occupied by the 
two extra electrons contributed by the arsenic atom that 
has been substituted for a gallium atom, which cannot 
be accommodated in the bonding states of the valence 
band. 

We conclude that the charge transition levels of the 
Vca and the As^-As should remain fixed relative to the 
valence band edge, while the donor levels of Asca, which 
possess a conduction band character, should be shifted 
up together with the conduction band states. In Fig. [5J 
we show the charge transition levels for these defects cor- 
rected by the above procedure, using the room temper- 
ature gap of 1.4 eV. The transition levels of Asc a are 
shifted to 1.0 eV and 1.1 eV, in fortuitously good agree- 
ment with the MCDA results identified with this defect 
in LT GaAs ; 14 ' 49 although the transitions are both about 
0.4 eV higher than those associated with Asc a in melt- 
grown GaAsi^S 

IV. SUMMARY 

In this work we have shown that larger supercells and 
a better k-space sampling than have been used in a 
large number of previous DFT defect calculations are 
required to give accurate results for the formation en- 
ergies, charge transition levels, defect state properties, 
and atomic structure and relaxation for the arsenic split 
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interstitial in GaAs. In particular, we find that 217-atom 
supercells are necessary to get good results for the charge 
transition levels and the dispersion of the deep defect 
state, even though the arrangement of atoms in the struc- 
ture is well converged in a 65-atom cell, particularly if one 
uses the finely spaced 2 3 MP k-space mesh. 

We have calculated formation energies for As split in- 
terstitials, Ga vacancies, and As antisites in As-rich GaAs 
using the larger supercells and better k-point sampling 
which we have determined to be necessary. Using these 
results, we find that the equilibrium concentrations of ar- 
senic interstitials will be substantially lower than equilib- 
rium concentrations of arsenic antisites in As-rich, n-type 



or semi-insulating GaAs. 
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